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ABSTRACT 

It has recently been suggested that, in the field, ~ 56% of Sun-like stars 
(0.8 M 0 < M* < 1.2 M 0 ) are single. We argue here that this suggestion may be in¬ 
correct, since it appears to be based on the multiplicity frequency of systems with 
Sun-like primaries, and therefore takes no account of Sun-like stars that are secondary 
(or higher-order) components in multiple systems. When these components are in¬ 
cluded in the reckoning, it seems likely that only ~ 46% of Sun-like stars are single. 
This estimate is based on a model in which the system mass function has the form pro¬ 
posed by Chabrier, with a power-law Salpeter extension to high masses; there is a flat 
distribution of mass ratios; and the probability that a system of mass M is a binary 
is 0.50 + 0.46 log 10 (M/M o ) for 0.08 M 0 < M < 12.5 M 0 , 0 for M < 0.08 M 0 , and 
1 for M > 12.5 M 0 . The constants in this last relation are chosen so that the model 
also reproduces the observed variation of multiplicity frequency with primary mass. 
However, the more qualitative conclusion, that a minority of Sun-like stars are single, 
holds up for virtually all reasonable values of the model parameters. Parenthetically, 
it is still likely that the majority of all stars in the field are single, but that is because 
most M Dwarfs probably are single. 
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1 INTRODUCTION 


The multiplicity statistics of Sun-like stars in the field have 
recently been re-evaluated by Raghavan et al. (20101, us¬ 
ing a volume-limited sample of 454 stars within 25 pc. They 
estimate that, for systems having a Sun-like primary, the 
multiplicity frequency is m s ( M 1 =M 0 ) = 0.44 ± 0.02rl 

The multiplicity frequency for systems having a primary 
of mass M x is defined as 


Wr) = 


B + T + Q+... 
S + B + T + Q + . 


( 1 . 1 ) 


where S is the number of single stars of mass M x , B is the 
number of binary systems having a primary of mass M x , 
T is the number of triple systems having a primary of mass 
M x , Q is the number of quadruple systems having a primary 
of mass M x , and so on. Implicit in the above statement is 
the fact that in reality one must consider a finite interval of 
mass, in order to have meaningful statistics. 


* E-mail: ant@astro.cf.ac.uk 

1 Here the subscript S records the fact that m s (M = M 0 ) is a 
system-property. “(M 1 =M 0 )” is to emphasise that — for practi¬ 
cal observational reasons - the convention is to estimate m g as a 
function of primary mass, M x , rather than system mass. We use 
m g for multiplicity frequency, and not mf as is normal, because 
mathematical expressions involving a double symbol like mf are 
inevitably confusing. 


Thus, the Raghavan et al. (20101 estimate of m s (M 0 ) 
indicates that the number of single Sun-like stars exceeds, 
by a factor of 1.28 ± 0.08, the number of Sun-like stars that 
are primaries in multiple systems. This is not quite the same 
as the inference made by Raghavan et al. (2010), “the ma¬ 
jority ... of solar-type stars are single”, and reiterated by 


Duchene & Kraus (20131 “a slight majority of all field solar- 


type stars are actually single”, since it takes no account of 
Sun-like stars that are secondaries, tertiaries, etc. in multiple 
systems, only those that are primaries. 


In the sequel we show that the majority of Sun-like 
stars may, in fact, be in multiples. To keep the analysis sim¬ 
ple we consider only single and binary systems, so we are 
concerned with evaluating the number of Sun-like stars that 
are secondaries in binary systems. Consideration of higher- 
order multiple systems would only strengthen our conclu¬ 
sion by bringing into the reckoning additional components 
that might be Sun-like. Section [2] presents the model of bi¬ 
nary statistics that we use, the ranges in which we allow the 
model parameters to vary, and the basic analysis. Section [3] 
presents the results, i.e. the fractions of Sun-like stars that 
are single, primary or secondary, and how these fractions de¬ 
pend on the parameter choices made. Section [4] summarises 
our conclusions. 
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2 MODEL PARAMETERS AND ANALYSIS 


2.2 The fraction of systems that are binaries 


In order to model the binary statistics of Sun-like stars and 
estimate what fraction are secondaries in binaries, we need 
to specify the distribution of system masses, the fraction of 
systems that are binaries (as a function of system mass), 
and the distribution of mass ratios in binary systems. 

We define the system mass to be M, and the masses 
of stars to be Mi, where i = 0 corresponds to a single star 
(hereafter a single), i = 1 corresponds to the primary in a 
binary (hereafter a primary), and i = 2 to the secondary in 
a binary (hereafter a secondary). In addition, we introduce 
the corresponding logarithmic variables 


a = lo §i n ^ ^ ’ (systems); 


Since 


M = M 0 KT = M 0 e^ , 
where t = ln( 10 ), we have 
dM 


(stars). 


= £M. 


dp 

The mass ratio of a binary system is q = M 2 /M 1 . 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 


2.1 The distribution of system masses 

We assume that low- and intermediate-mass systems (p < 
p c , see below) have a log-normal mass distribution, as pro¬ 
posed by Chabrier (20051, 

djV 


(2tt) 1 /2 c 


exp 


(M Ms ) 

2a| 


(2.5) 


where the subscript ln is for log-normal. Chabrier (2005) 


estimates that p s = —0.60 and a s = 0.55; the normalisa¬ 
tion in [Chabrier (2005) is different because he uses physical 
units. We assume that these values are accurate to ±0.05, 
i.e. —0.65 < p s < —0.55 and 0.50 < a s < 0.60. 

For higher-mass systems, M > M c , we adopt a power- 
law distribution of masses, with exponent a, i.e. 


dJV 
dM , 


K 

JuZ 


M_ 


( 2 . 6 ) 


0 \ ©V 

where the subscript PL is for power-law. It follows that 

dJV 

^ i_ 

The default exponent is a = 2.35, as first estimated by 


= K exp{—(a — 1 )£p} . 


(2.7) 


Salpeter (1955), and we assume this is accurate to ±0.35, 


i.e. 2.00<a<2.70. 

We require that the join between the two distributions 


be smooth, i.e. at p c the distributions (Eqns. 2.5 & 2.71 and 
their slopes should be equal, so 


Me = Ms + ( Q - l )tc 


K = 


(2t r)!/2 c 


exp< (a — l)lp c — 


(Me Ms) 

2°s 


( 2 . 8 ) 

(2.9) 


With the default values of the model parameters, the join 
occurs at p c ~ 0.34, corresponding to a system mass of M c ~ 
2.2 M„. 


Since the distribution of mass ratios appears to vary slowly - 
if at all with primary mass in the range 0.2 M 0 <M X <2 M 0 


__ _______ _ '' ^ _Li 

(Janson et al.|2012 Raghavan et al.|2010 Rcggiani & Meyer 


20131, the fraction of systems that are binaries must depend 


on system mass in a similar way to the dependence of multi¬ 
plicity frequency on primary mass, but displaced to slightly 
higher masses. We therefore assume that the fraction of sys¬ 
tems that are binaries is given by 

f 0 , n<~Po/Pi\ 

Pin) = { Po+PiV, -Po/Pi <M<(1 - Po)/Pi', (2.10) 

l 1, m> ( 1 -&)/&. 

The justification for adopting this functional form is given 
in §2.6| where we also derive the default values of (3 0 and /3j, 
and their ranges, by fitting the observed run of multiplicity 
frequency, m s against primary log-mass, p, (see Fig. [l]). 


2.3 The distribution of mass ratios 

Estimates of mass ratios are sufficiently uncertain, and the 
samples subject to such severe selection effects, that the ap¬ 
propriate distribution is poorly constrained. A convention 
has emerged (e.g. Reggiani & Meyer 2013 Duchene & Kraus 
20131 that - at least in the first instance - the distribution 
should be fitted with a power-law (exponent 7 ), possibly 
with a minimum mass ratio (g MIN ). We therefore adopt the 
distribution 


dP 

dq 


0, 


q<q* 


(7 +1) (1 - 4 ™ 1) ) - q > 9 iv 


( 2 . 11 ) 


The default values are q Mm = 0.00 and 7 = 0.00, i.e. a flat 
distribution with no minimum (e.g. |Raghavan et al.]|2010| 
Janson et al.|20l2 |. 

Finite q Mm reflects the possibility that Sun-like pri¬ 
maries may eschew brown-dwarf or very low-mass secon¬ 
daries (the so-called Brown Dwarf Desert). We therefore 
consider the range 0.00 < <j MIN < 0 . 10 . 

Regarding the range of 7 , we adopt —0.60 < 7 < 2.00 
(see Duchene fe Kraus|2013 and references therein). 

Negative 7 means a preference for low-mass compan¬ 
ions, and 7 = —0.60 means that (with g MIN = 0) ~ 50% 
of systems have q < 0.2. This might be appropriate for 
intermediate-mass and/or wide systems. However, the ob¬ 
servational evidence for a predominance of such low q values 
may actually refer to low-mass tertiary components orbiting 
unresolved binaries. 

Positive 7 means a preference for a companion of com¬ 


parable mass (or “twin”, e.g. Lucy 2006), and 7 = 2.00 means 
that (with q Mm = 0) ~ 50% of systems have q > 0.8. How¬ 
ever, the observational evidence for a predominance of such 
high q values may reflect selection effects. It is mainly seen 
in systems at the extremes of primary mass, where the sta¬ 
tistical limitations of the observational data are most severe, 
and our concern here is not with such extreme masses. 
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Source 

m s (M 0 ) 

Pn 

Pi 

MMJ 

/*(M 0 ) 

jRaghavan et al. (j2010|) 

0.44 

0.50 

0.46 

0.535 

0.461 

Duquennoy & Mayor) Jl991|) 

0.58 

0.64 

0.52 

0.667 

0.450 


Table 1. Column 1 gives the source of the estimate of the multiplicity frequency for systems with Sun-like primaries, and 

Column 2 gives the value of m s (M 0 ). Columns 3 and 4 give the values of the parameters /3 0 and /3 1 in the expression for the probability 
that a system of given log-mass is a binary (Eqn. 2.10). Columns 5 and 6 give the fraction of Sun-like stars that are in binaries, b* (M 0 ), 
and the ratio of Sun-like secondaries to Sun-like primaries, /* (M 0 ). 


2.4 The distributions of singles, primaries and 
secondaries 


With these definitions, the mass distribution of singles is 


dAf 


dUo 


dAf 


dfi 


~t~ = -7- U-/3(mo)} ; 


( 2 . 12 ) 


the mass distribution of primaries is 


dAf 

dUi 


i+l°g 10 (2) 


/ 


dV . 


i 


dP 

dq 


t{q+ 1) 


dn\ 

q = io ( e-ei)_i 

(2.13) 


Rizutto et al Mason et al 



fj, l = log (Mj/Mg) (primary log-mass) 


and the mass distribution of secondaries is 

/j,=oo 


dAf 
dfi 2 


/ 


M=M 2 + lo g 10 ( 2 ) 


dV 

dfi 


Pin) 




q=('l0 ( ' i -' i 2)-iy 


dfi. 


(2.14) 


The last terms in the integrands of Eqns. (2.131 and (2.141 
represent, respectively, \dP/dfi\ lll and |<9P/<9mIm 2 - 


2.5 Statistics 


To obtain the binary statistics of stars having mass A/*, 
we set Mo = Hi = /+ = l°g lo (Af*/M 0 ) in Eqns. (2.12|, 
(2.131 and (2.14|, and evaluate the integrals numerically. 
The multiplicity fraction for systems having primaries of 
this mass is then 


m s (A/*) = 


dV /dfi 1 


dV /dfi 0 + dAf /dfi 1 
but the fraction of such stars that are in binaries is 
dAf/dfi 1 + dAf /dfi 2 


6* (A/*) = 


dAf/dfi 0 + dAf /dfi 1 + dAf /dfi 2 


(2.15) 


(2.16) 


We can also compute the ratio of stars of mass Ad* that are 
secondaries, to stars of mass Af* that are primaries, 


f*(MA 


dAf /dfi 2 
dAf /dfi 1 


/ (b* -m s ) \ 

V m s (i- 6 *)/ ’ 


(2.17) 


Figure 1. The variation of multiplicity frequency, m s , with pri¬ 
mary log-mass, /x 1 = log 1Q (M x /M 0 ). The hatched acceptance 
boxes represent observational estimates fro m|Close et al7| (|20Q3|) , 
|Basri fc Reiners] (|2006|), |Fischer fc Marcy| <|1992[), |Janson et al.| 
(|2012|) , |Duquennoy fc Mayor | (| 1991 1), |Kouwenhoven et al.| <|2007[), 

|Rizzuto et al.| < |2013[ ), |Preibisch et al.| ( |1999[ ), and | Mason et al.| 

( 1998J) ; the width of each box represents the approximate range 
of masses considered, and the height of each box represents the 
range of m s deduced. Note that, for the higher mass ranges (i.e. 
Kouwenhoven et al.||2007| |Rizzuto et al.|[2013[ |Preibisch et al.~] 
1999||Mason et al.|1998[ ), there are only lower limits on ra s . The 
lower (bolder) line is the best fit when we require the model to go 
through the |Raghavan et alT] ( |2Q10| estimate of m s for systems 
with Sun-like primaries; this fit defines the default values of the 
parameters (3 0 and (3 1 (see Table [2j). The upper (feinter) line is the 
best fit when we require the model to go through the |Duquennoy | 
|fe Mayor] < |1991[ ) estimate of m s for systems with Sun-like pri¬ 
maries; evidently this line fits the other data more comfortably. 


2.6 Constraining /3 0 and /3 X 

We identify the default values of /3 0 and p 1 by requiring 
that the model (i) reproduce accurately the multiplicity fre¬ 
quency of Sun-like stars inferred from observation by |RaghaT] 
van et al. (2010), and (ii) reproduce as closely as possible 


the run of multiplicity frequency with primary mass derived 


from other observational studies by Close et al. (2003), Basri 


& Reiners (20061, Fischer & Marcy ( 

1992), Janson et al. 

(2012 

1, Duquennoy & Mayor (1991), 

Kouwenhoven et al. 


(2007|), Rizzuto et al. (2013), Preibisch et al. (1999) and Ma¬ 


son et al. (1998). The resulting default values are /3 0 =0.50 


and /3 X =0.46, so the probability that a system of mass M 
is a binary becomes 


this is the fractional increase in the number of stars of mass 
M* in binaries that derives from taking account of secon¬ 
daries. The subscripts ★ on b and / record that these are 
stellar properties. 


P{M) 


0, M < 0.08 M 0 ; 

0.50 +0.46 log lQ (A//M 0 ), 

0.08 M 0 <M < 12.5 M 0 ; 
1, M > 12.5 M 0 . 
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MODEL 

MODEL 

DEFAULT 

RANGE 

dm g 

db * 

dU 

FUNCTION 

PARAMETER, X 

VALUE, X Q 

AX 

dX 

dX 

dX 

DISTRIBUTION 

Ms 

-0.60 

±0.05 

0.13 

0.18 

0.27 

OF SYSTEM 

°S 

+0.55 

±0.05 

0.32 

0.44 

0.75 

LOG-MASSES 

a 

+2.35 

±0.35 

0.00 

-0.02 

-0.14 

BINARY 

ho 

+0.50 

±0.05 

0.90 

0.90 

-0.20 

FRACTION 

hi 

+0.46 

±0.05 

0.07 

0.09 

0.15 

DISTRIBUTION 

^MIN 

+0.00 

+0.10, -0.00 

-0.07 

0.03 

0.59 

OF MASS RATIOS 

7 

+0.00 

+2.00, -0.60 

-0.03 

0.01 

0.26 


Table 2. Column 1 gives the model parameters. Columns 2 through 4 give their symbols, their default values and their ranges. Columns 
5 through 7 give the derivatives of, respectively, the multiplicity frequency of Sun-like stars, m s (M 0 ), the fraction of Sun-like stars that 
are in binaries, 6*(M 0 ), and the ratio of Sun-like secondaries to Sun-like primaries, /*(M 0 ), with respect to each model parameter in 
turn; all the derivatives are evaluated for the default values of the model parameters. 


The corresponding fit to the observational estimates is illus¬ 
trated by the lower (bolder) line on Fig. [I] A brief discussion 
of this plot is appropriate. 

The model equations and default parameters describing 
both the distribution of system masses (i.e. Chabrier log¬ 
normal at low and intermediate masses, with /x s = —0.60 
and u s =0.55, plus Salpeter power law, with negative slope 
a = 2.35, at high masses), and the distribution of stellar 
mass ratios (flat, 7 = 0 . 00 , with no minimum, q M1K = 0 . 00 ) 
appear to be the natural default choices; we will explore the 
consequences of varying the model parameters, both individ¬ 
ually and collectively, in f|3] The justification for adopting 
Eqn. (2.101 for the fraction of systems that are binaries as a 
function of system mass, and hence the choice of the model 
parameters /3 0 and hi , is more adhoc. 

From the observational data presented in Fig. [T] it ap¬ 
pears that, for — m s increases approximately 

linearly with fi 1 , and hence a linear relation between /? and 
fj, (i.e. Eqn. 2.10 1 seems the simplest option to explore. How¬ 
ever, it is also clear from Fig. |l| t hat Eqn. (2.101 gives a 
much better fit when the Duquennoy & Mayor (1991) ac¬ 
ceptance box is invoked for systems with Sun-like primaries 
(upper, feinter line) than whan the Raghavan et al. (2010) 
acceptance box is invoked (lower, bolder line). Even if the 
Fischer & Marcy (|1992 1 and Duquennoy fe Mayorj (19911 ac¬ 


ceptance boxes are discounted, the Raghavan et al. (20101 
acceptance box is hard to reconcile with the others with¬ 
out introducing a more complex function for which 

the slope has a distinct minimum around /+ ~ 1. Since - as 
far as we are aware - no physical explanantion for such a 
minimum has been advanced, we avoid this complication. It 
would certainly be convenient if, when in future the multi¬ 
plicity frequencies for systems with non-Sun-like primaries 
are evaluated more accurately, the values are reduced (say 
by ~ 20%) to bring them into line with Raghavan et al. 


(20101, or if the Raghavan et al. (20101 estimate is revised 
upwards. 

In this context, it is important to note that Sun-like 
singles derive from systems with /.t = 0, Sun-like primaries 
from systems with 0 </j,< log 1Q (2), and Sun-like secondaries 
from systems with log iQ (2) < /x < 00 . In addition, with any 
sensible choice of the model parameters, the system mass 
function is falling quite rapdly with increasing /x for /x > 0 . 
Thus the important range of applicability of Eqn. (2.101 is 
0</x<0.6. 

At low /x (/x < —ho/hu M < 0.08M 0 ), /3 has to be 
set to zero to avoid non-physical predictions, and this has 


the consequence that the multiplicity frequency falls to zero 
for jU x < —1.3 (i.e. M 1 <O.O5M 0 ). The observed multiplicity 
frequency appears to be very low for such low-mass primaries 
(only one of the systems reported by Close et al. (2003) 
has M 1 <O.O5M 0 ), but it is probably not zero. We explore 
the effect of adopting model parameters that increase the 
multiplicity frequency of brown dwarfs and very low-mass 
stars in © However, it is important to note that the value 
of (3 at these low masses is irrelevant to the statistics of Sun¬ 
like stars, because a low-mass core, M < M 0 , cannot spawn 
a Sun-like star. 

At high n (fi > (1 — ho)/hi> M > 12.5M 0 ), /3 has to 
be set to one to avoid non-physical consequences. In fact the 
acceptance boxes due tojKouwenhoven et al.|(|2007|),|Rizzuto| 


et al. (2013), Preibisch et al. (1999) and Mason et al. (1998 1 


are all lower limits. At these high primary masses, the mul¬ 
tiplicity frequency is an inadequate measure of multiplicity, 
because higher-order multiple systems become increasingly 
important at high primary mass, and the multiplicity fre¬ 
quency does not distinguish between higher-order multiples 
and binaries (fHubber & Whitworthj20051. The pairing fac¬ 
tor (number of orbits per system) or companion frequency 
(mean number of companions) would be more appropriate 
measures. In addition, the notion of liigh-mass field stars 
is fraught, because the highest-mass stars do not live long 
enough for their birth clusters to disolve completely into the 
held, and therefore many high-mass stars in the held are 
runaways, which have been ejected in violent A-body in¬ 
teractions, and therefore seldom, if ever, have companions. 
The multiplicity statistics for high-mass stars, which point 
to very high pairing factors and companion frequencies, tend 
actually to pertain to stars that are still intimately involved 
with their birth clusters. 


3 RESULTS 

3.1 The default solution 

For Sun-like stars and the default parameter set, we obtain 
m s (M 0 ) = 0.440, b*(M 0 ) = 0.535 and /*(M 0 ) = 0.42. In 
other words, although for every ~ 56 single Sun-like stars, 
there are only ~ 44 Sun-like primaries in binary systems, 
there are also ~ 20 Sun-like secondaries in binary systems, 
hence a total of ~66 Sun-like stars in binaries. 
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3.2 The effect of varying the model parameters 


Figs. E and [3] represent the (m s (M g ), 6 *(M 0 ))-plane in 
the vicinity of the solution obtained with the default pa¬ 
rameters (hereafter the default solution). The black dot at 


m s (M 0 ) = 0.440 and 6 *(M 0 ) = 0.535 on Fig. [2 
the default solution. The vertical dashed lines 

represents 
emark the 

limits on m s (M 0 ) obtained by 
0.42 < m s (M 0 ) < 0.46. The 

Raghavan et al. 

| 2010 ), viz. 

lorizontal dotted line is at 


&*(M 0 ) = 0.50; for solutions below (above) this line more 
(less) than 50% of Sun-like stars are single. 

On Fig. [2} the lines passing through the default solution 
show how the solution changes if one of the model parame¬ 
ters is varied. The arrow at one end of a line indicates the 
direction in which the parameter increases, and the symbol 
for the parameter in question is given beside this arrow. The 
lines for p , s , cr s and /% lie almost on top of one another; each 
of these lines extends approximately the same distance on 
either side of the default solution, and therefore the length 
can be inferred from the position of the corresponding arrow 
(the line for a s is the longest, and that for /% the shortest). 


Over the range of system log-masses that contributes 
Sun-like stars, 0 < p, < 00 , the mass function is decreasing 
with increasing fi, i.e. dAT/dfi <0. Increasing p s and/or in¬ 
creasing (T s and/or decreasing a increases dAf/dp < 0, i.e. 
reduces the downward slope of the system mass function, 
so there are then more Sun-like primaries relative to singles 
(increased m s (M 0 )) and more Sun-like secondaries relative 
to primaries (increased 6 *(M 0 )). Conversely, decreasing p s 
and/or decreasing <r s and/or increasing a decreases m s (M 0 ) 
and 6 *(M 0 ). Increasing a from its minimum value (a = 2.00) 
has an ever diminishing effect, because the switch from the 
log-normal system mass distribution to the power-law dis¬ 
tribution occurs at ever increasing p c (see Eqn. 2.8 1 ; conse¬ 
quently the line for increasing a does not extend far beyond 
the default solution on Fig. [2} once a >2.35. 


Increasing /3 0 and/or /3 1 above the default values in¬ 
creases the number of primaries and secondaries, relative to 
singles, and hence increases both m s (M 0 ) and 6 *(M 0 ); in¬ 
creasing f3 0 has a larger effect than increasing /3 1 because the 
binary statistics of Sun-like stars are dominated by systems 
with 0 < /r<0.6, i.e. relatively small p. Conversely, decreas¬ 
ing /3 0 and/or /3 1 below their default values decreases both 
m s (M 0 ) and b*(M 0 ). 


Increasing q MIN and/or 7 above the default values 
squeezes the distribution of systems involving Sun-like pri¬ 
maries towards the upper end of the available range (0 < 
p < log l 0 ( 2 )), i.e. down the falling system log-mass distri¬ 
bution, thereby reducing the number of Sun-like primaries. 
At the same time, it squeezes the distribution of systems 
involving Sun-like secondaries towards the lower end of the 
available range (log 10 ( 2 ) < p < 00), i.e. up the falling sys¬ 
tem log-mass distribution, thereby augmenting the number 
of Sun-like secondaries. Consequently m s (M 0 ) is reduced, 
but 6 *(M 0 ) is augmented. q MIN cannot be reduced below 
the default value, but decreasing 7 below the default value 
augments m s (M 0 ) and reduces 6 *(M 0 ). 

Table [2] gives the parameter symbols; their default 
values and prescribed ranges; and the partial deriva¬ 
tives, dm s (M Q )/dX, i%*( M @ )/dX and <9/*(M 0 )/<9A, 
where X is one of the model parameters (i.e. X = 


p s , a s , (3 01 q MIN , 7)5 an< i the derivatives are evalu¬ 
ated for the default parameter set. 

Fig-! displays the distribution of solutions when the 
model parameters are varied simultaneously and randomly. 
This plot is produced by generating ~ 1.8 x 10 9 different so¬ 
lutions (~ 10 9 of which fall within Fig. [ 3 ]) . For each solution, 
the value of parameter A' is chosen by generating a random 
Gaussian deviate, Q, from a distribution with mean 0 and 
standard deviation 1, and then putting X = A' 0 +Q AX, 
where the values of X Q and AA' for each model parameter 
are given in Table [ 2 ] (Columns 3 and 4); with this procedure, 
~ 32% of parameter values fall outside the range A 0 ± AA. 

If we consider all the solutions generated in this way 
(including those that fall outside Fig. [3|, ~ 21% of them 
give /* <0.50. However, the majority of these solutions in¬ 
volve very low values of /3 0 , and therefore they also deliver 
m s (M 0 ) < 0.42 (i.e. below the range estimated by 
[et al.pTlO). 


Raghavan 


If we limit consideration to solutions that satisfy the 
constraints calculated by Raghavan et al. (2010), i.e. 0.42 < 
m s (M 0 ) < 0.46 (including those that fall outside Fig. |3j) , 
then ~ 3% of the allowed solutions give /* < 0.50, because 
these solutions require both very low /3 0 and very low 7 . 
We conclude that it is rather unlikely that the majority of 
Sun-like stars are single. 


3.3 Other masses 

If we accept the default parameters as a reasonable represen¬ 
tation of the binary statistics of low- and intermediate-mass 
stars, we can estimate the fraction of such stars that are 
single. For example, for M* < 0.70 M 0 , > 50% of stars are 
single. Thus most M Dwarfs are single. 


4 CONCLUSIONS 


We have developed a model to estimate the fraction of Sun¬ 
like stars that are secondaries in binary systems. Our model 
(which only considers singles and binary systems, and there¬ 
fore gives a lower limit to the fraction of Sun-like stars that 
are in multiples) invokes a log-normal distribution of sys¬ 


tem masses (Chabrier 


2005) with a power-law tail at high 


masses (|Salpeter[ 19551; a flat distribution of mass ratios 


(e.g. |Raghavan et al.||2010{ |Janson et al.|[2012| |Reggiani fe| 


Meyer 20131; and the probability that a system with mass 
M is a binary, 


13 = 0.50 + 0.47 log (M/M 0 ). 


We find that, even if the multiplicity frequency is as low 


as the recent estimate of Raghavan et al. (2010), 0.42 < 


m s (M 0 ) < 0.46, and even if we vary the model parame¬ 
ters significantly from their default values, when account is 
taken of secondaries (and higher-order components in mu- 
tiple systems), the majority of Sun-like stars are probably 
in multiple systems. Our model predicts that the fraction of 
stars that are in multiple systems is a monotonically increas¬ 
ing function of stellar mass, and that the majority of stars 
with mass M* >O.7M 0 are in multiple systems. Conversely, 
most M Dwarfs, and hence most stars overall, are single. 
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Figure 2. The interdependence of the multiplicity frequency for 
systems having Sun-like primaries (m s (M 0 ); abscissa) and the 
fraction of Sun-like stars that are in binary systems (5*(M 0 ); 
ordinate). The filled circle near the centre represents the combi¬ 
nation obtained with the default parameters (see Table [2). The 
lines track the changes in m s (M 0 ) and 6*(M 0 ) that occur when 
the model parameters are changed individually between the limits 
in Table [ 2 ] (for example, when fi s varies from -0.65 to -0.55). In 
each case, an arrow at one end of the line indicates the direction in 
which the parameter increases, and the symbol for the parameter 
in question is always located beside this arrow. The lines for fi s , 
o' s and /3 X are almost exactly on top of one another. Since each 
of these lines extends almost the same distance on either side of 
the default point, their lengths can be estimated from the loca¬ 
tions of the corresponding arrows; the longest line is for <j s and 
the shortest for f3 ± . The two vertical dashed lines mark the range 
of m s (M 0 ) deduced by Raghavan et al. ( 2010). The horizontal 
dotted line is at 6*(M 0 ) = 0.500, so all solutions abve this line 
represent situations in which the majority of Sun-like stars are in 
multiple systems. We see that no individual parameter can reduce 
5* (M 0 ) below this line without also taking m s (M 0 ) outside the 
range derived by Raghavan et al. (2 010[ >. 


Figure 3. The likelihood of different combinations of m s (M 0 ) 
(the multiplicity frequency for systems having Sun-like primaries; 
abscissa) and 5*(M 0 ) (the fraction of Sun-like stars that are in 
binary systems; ordinate). Different random sets of model param¬ 
eters have been generated by picking a value for each parameter 
from a Gaussian distribution having mean equal to the default 
value and standard deviation equal to the specified range; thus, 
for any parameter there is a ~ 32% chance that it falls outside the 
range X Q ± AX. The grey-scale and contours delineate the rela¬ 
tive probability of different combinations of m s (M 0 ) and 6*(M 0 ) 
that fall within the bounds of the figure. The contours are at 15%, 
30%, ... 90% of the peak probability. If we consider all the combi¬ 
nations generated (including those that fall outside the bounds of 
the figure), there is a ~21% chance that the majority of Sun-like 
stars are single. However, this result is dominated by parame¬ 
ter sets with very low /3 X which also deliver m s (M 0 ) well below 
the minimum estimated by |Raghavan et alT] ( |2010| ), i.e. the ex¬ 
tension into the bottom lefthand corner of the figure. If we limit 
our consideration to the sets that deliver ra s (M Q ) in the range 
deduced by |Raghavan et al. | ( j20io| , the chance that the majority 
of Sun-like stars are single falls to ~3%. 
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